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Abstract 

^ ' The autocorrelation function of the force acting on a slow classical sys- 

>D , tern, resulting from interaction with a fast quantum system is calculated 

following Berry-Robbins and Jarzynski within the leading order correction 
to the adiabatic approximation. The time integral of the autocorrelation 
function is proportional to the rate of dissipation. The fast quantum sys- 
tem is assumed to be chaotic in the classical limit for each configuration of 
^v [ the slow system. An analytic formula is obtained for the finite time inte- 

gral of the correlation function, in the framework of random matrix theory 
(RMT) , for a specific dependence on the adiabatically varying parameter. 
Extension to a wider class of RMT models is discussed. For the Gaussian 
unitary and symplectic ensembles for long times the time integral of the 
Q ' correlation function vanishes or falls off as a Gaussian with a character- 

Cu , istic time that is proportional to the Heisenberg time, depending on the 

1-^ ' details of the model. The fall off is inversely proportional to time for the 

Gaussian orthogonal ensemble. The correlation function is found to be 
dominated by the nearest neighbor level spacings. It was calculated for 
1 J a variety of nearest neighbor level spacing distributions, including ones 

rS I that do not originate from RMT ensembles. The various approximate 

^ . formulas obtained are tested numerically in RMT. The results shed light 

on the quantum to classical crossover for chaotic systems. The implica- 
tions on the possibility to experimentally observe deterministic friction 
are discussed. 

PACS: 05.45.Mt, 05.45.Ac, 05.45.-a 

1 Introduction 

Dissipation of energy from a physical system to a thermal bath takes place as a 
result of a fluctuating force that acts on the system because of its coupling to the 
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bath. The dissipative friction force is proportional to the correlation function 
of the fluctuating force. The force resulting from coupling to a chaotic system 
is also fluctuating. The question that will be studied in the present paper is on 
what time scales it leads to friction and what is the relation of this friction to the 
autocorrelation function of the fluctuating force. An example is slow particle 
coupled to a fast particle so that the motion of the fast particle is chaotic for 
each position of the slow particle. 

Various models for dissipation of energy from a slow particle by a fast one 
have been developed. To our knowledge the first models of this type were intro- 
duced in the context of nuclear physics Q . In particular, a model where particles 
move within a region bounded by a deforming boundary, modeling the nuclear 
surface was studied [^ . The energy transferred between the boundary and par- 
ticles enclosed inside was calculated classically and quantum mechanically in 
the framework of some approximations. Recently some detailed numerical sim- 
ulations were performed along these lines and the regime of validity of various 
approximations was tested 0] . The dissipation for a wide class of model systems 
was explored by Wilkinson and Austin [4-6] in the framework of random matrix 
theory (RMT), relevant for a situation where the motion of the fast particles 
is chaotic. The relation to Landau-Zener tunneling was also studied in these 
works. A different RMT model was studied by Mizutori and Aberg m. In all 
these studies dissipation was found. In addition, a different approach, aimed at 
emphasizing the relation to many body problems was introduced Isl. A unified 
picture of many of these models has recently been presented by Cohen M . A 
systematic investigation of the interaction of a slow system with a fast one is 
possible with the help of multiple scale analysis. Under such conditions Ott 
demonstrated fI3| that the phase space volume enclosed by the energy surface 
of the fast particle is an adiabatic invariant, namely its change is much slower 
than that of the fast particle Hamiltonian. It has been demonstrated for various 
conditions that it is indeed an adiabatic invariant [|ll| . In the present paper we 
study the behavior of a slow particle that is coupled to a fast chaotic system. A 
model for such a system, that is quite general, and has been studied by Berry 
and Robbins (BR) iQ and by Jarzynski [^ in the framework of multiple scale 
analysis, is defined by the Hamiltonian: 

The phase space coordinates of the slow particle are (P, R) and its mass is M. 
For simplicity it is coupled only through its position to the fast system whose 
phase space coordinates are z = (p,r). The latter system has the property 
that if R is kept fixed it is fully chaotic. The crucial feature of the system we 
wish to study in this work is that it exhibits a wide separation of time scales- 
the evolution of the fast system, characterized by the time scale Tfag+i is so 
rapid that it explores all of the phase space accessible to it energetically before 
the slow particle, characterized by the time scale Tgi„ , moves appreciably. 



The adiabaticity parameter is £ ~ -^fast/^slow One way to realize this is 
to couple two particles with a mass ratio of ra/M — e^ ^ 1, as one can see 
by rescaling the equations of motion. We now turn to analyze the dynamics 
generated by the Hamiltonian dl^) with the approximation that the slow particle 
evolves under the influence of the average force exerted on it by the fast system, 
which can be treated as a system described by a slowly varying Hamiltonian. 
The time dependence of this Hamiltonian is determined by the dynamics of the 
slow particle. First the classical dynamics is outlined and later the quantum 
mechanical behavior is summarized. In the case of (yj) the average force is given 
by: 

F(Ta) = - fdzp{z,Ta)dRh{z,Tl{Ta)), (2) 

where p(z,,Ta) is a normalized probability density in the fast particle phase 
space. A formalism that includes fluctuations was developed by Jarzynski p3| . 
The results of this paper do not depend on these fluctuations and therefore the 
formalism of BR will be used. The probability density satisfies the Liouville 
equation: 

e^P(z,Ta) = ^h{z,'R{Ta)),p{z,Ta) |^, (3) 

written in a way that emphasizes that the evolution of the fast system is indeed 
on a much shorter time scale than the time scale on which the fast Hamiltonian 
changes. {}z denotes Poisson brackets taken with respect to z. With the aid of 
the multiple scale expansion: 

oo 

P(z,T-a) =^eVi(z,Ta), (4) 

1=0 

Berry and Robbins |1^ were able to calculate the force acting on the slow 
particle up to first order in e: 

F«Fo + £Fi. (5) 

To leading order, the force is given by the classical analogue of the Born- 
Oppcnheimcr force: 

Fo,iTa)^~dR,E{-R), (6) 

where E (R) is chosen such that the phase space volume enclosed by the energy 
surface of the fast particle, fl {E (R) , R), is constant. The leading correction to 
Fo includes a velocity dependent force: 



Fuira) = - > K,,R,; K,, = S^^a^ 



J:{E,R)-I,j{E,II) 



, (7) 

E=E{R} 



where !](£;, R) = dE^{E,R) and: 

/y {E, R) = 2 / di' Cy (^, R; t') ; 
Jo 

a,(S,R;t') = (afl,/i(zt,(z,R),R)afl,^./i(z,R) 



B.R 



In the last equation: 

/dz pQ{z,Ta)... = E-^(^,R) /dz 5(£;-/i(z,R)). 



/_E,R 



(8) 
(9) 

(10) 



denotes the microcanonical average and /i (z, R) = /i (z, R) — E{IV). FinaUy, 
Zt/(z,R) in (||) is the classical trajectory obtained by integrating the equations 
of motion generated by h (z, R) with fixed R backwards from z for the time t' . 
In addition to the velocity dependent force (M), and to the same order in £, there 
is a force that does not depend on velocity JMl. This force can be expressed 
as a gradient of a time dependent potential, and therefore is a correction to 
the Born-Oppenheimer force (0). Because of the form of Fi(Ta), it describes 
two qualitatively different forces. The first of these is geometric magnetism, 
and it is related to the antisymmetric part of Kij . This force has been studied 
analytically in the systems under discussion by BR | p2| and numerically by 
Berry and Sinclair fl^. The second force is associated with the symmetric 
part of Kij and is related to deterministic friction. This force has been studied 



in g, 16, |T2|, |l^. A central question that can be addressed at this point is 
under which conditions does the slow particle feel friction due to the velocity 
dependent force Fi. In order for this to happen Kij of (|^) has to have a positive 
definite symmetric part. 

The behavior when the fast system is quantum mechanical, to first order 
in e, has also been studied by BR ||l^. They found that in this case Kij 
is an antisymmetric tensor, meaning that the system exhibits only geometric 
magnetism and no friction. This difference is a result of the discreteness of 
the quantum spectrum. The quantum correlation function corresponding to the 



symmetric part of Cij of (g|) is: 






9R,h 



OrM 



n ) cos 



T {En — Ejn) 



(11) 



and the infinite time integral over it vanishes p2[ . In order to understand how 
the crossover between the classical and quantum behavior occurs, it is instructive 
to calculate the integral of the correlation function over a finite time. Since 
the discordance between the quantum and the classical models appears in the 



symmetric part of Kij it can be studied for the case where R is replaced be a 

scalar, time dependent parameter, X. In this case ^< Cij{n;t) + Cji{n;t) I will 

be denoted by C{t) for simplicity {Ri and Rj on the RHS of (O) will then be 
replaced by the scalar parameter X). Following BR we assume in the calculation 
that the initial state is an eigenstate of the Hamiltonian, |n). It was verified by 
BR that their result holds also if the initial state is a mixture. Our calculation 
can also be extended to a mixture leading to the same results. If the initial 
state is a pure state but not an eigenstate of the Hamiltonian one can check 
that within the assumptions of the paper the results are similar to the ones 
found if the initial state is an eigenstate of the Hamiltonian. The finite time 
integral that should be calculated then is: 

I{t) = f C{t')dt'. (12) 

Jo 

The correlation function and its integral may depend on the initial state n. This 
dependence has been suppressed in the notation for simplicity. In what follows, 
it will become clear that it is not important for the results of the present paper. 
Taking the classical limit ?i — > for any finite t and then the limit i — > cx) should 
result in a non- vanishing value of /(oo), while for any finite value of Ti, /(oo) 
should vanish. The friction on the time scale t is proportional to I{t) as can 
easily be inferred from (R) and (B) . The experimental meaning of this statement 
will be clarified in what follows. In order to understand the mechanism of this 
discordance, BR studied a model correlation function where the levels were 
equally spaced. They found that the function is periodic in time with period tp = 
h/AE, where AE is the level spacing. Moreover, in the classical limit, which 
in their model corresponds to taking tp —f oo, C{t) approaches the classical 
correlation function. For systems whose classical dynamics is chaotic, the energy 
levels are not equally spaced, but rather are distributed according to RMT IJITJ . 
The long time behavior of /(f) is determined by the levels nearest to n, namely 
nil, as can be seen from (|l l|) . Therefore, one may expect behavior different 
from the one found for the equally spaced spectrum. The natural question to 
ask is whether there is a characteristic time scale for the crossover between the 
quantum behavior of the integral I{t) and its classical behavior. The most naive 
answer to this question is that the characteristic time scale is the Heisenberg 
time, because it is the only time scale in the problem, and it is on this time scale 
that the quantum to classical crossover usually takes place. On the other hand, 
one can argue that there is no time scale for this crossover at all [Q . In RMT the 
probability for two consecutive levels to be separated by a distance s behaves like 
s^ for small spacings [19-22]. Consequently, {I pit)) ~ t~^ for long times, where 
here (. . .) denotes the RMT ensemble average, and Ip is the integral (O) for 
some (3. The answer given by the analysis presented in this paper is surprising. 
For the Gaussian orthogonal ensemble (GOE) {(3 = 1) one indeed finds that 
{Ip{t)) decays like 1/t, but for the Gaussian unitary ensemble (GUE) (/3 — 2) 



one finds that it decays like a Gaussian with a characteristic time proportional 
to the Heisenberg time or vanishes after the Heisenberg time depending of the 
parametric dependence on X. The integral //3(i) was also calculated for other 
values of (3. Why is the nature of the decay of I pit) important? There is the 
quantum-classical discordance that has already been mentioned, and one would 
like to analyze the scale that is required to observe the crossover between the 
regimes. It is relevant for some experiments, that will be mentioned below. 
In addition, there is the issue of the relevance of such an effect. The time 
over which the correlation function decays should be compared with other time 
scales present in the specific system studied. One such time scale is T2 ^ e~^, 
which is the time scale for the breakdown of the first order of the multiple scale 
analysis. Non-perturbative effects, such as Landau-Zener tunneling, become 
important on a time scale of T^^. In realistic experiments there is also the time 
scale for quantum decoherence T^j,. In order to observe the classical to quantum 
crossover discussed in the present work {Ip{t)) should exhibit substantial decay 
for t <C min {T2, Tlz) and of the order of T^. 

The model discussed in the present work is relevant for some experimental 
situations. Consider for example a molecular beam prepared in a classical con- 
figuration, where initially many levels are substantially populated. The beam 
travels in a slowly varying field p3[. Consequently the internal dynamics in 
the molecules is in a slowly varying potential. On short time scales the behav- 
ior is classical, the integral of the correlation function is positive and energy 
is absorbed in the motion of the internal degrees of freedom. On longer time 
scales the integral of the correlation function decays to zero and therefore one 
realizes that actually no energy is absorbed by the molecules. The outcome of 
the experiment depends on the time scale of decoherence, T^, that is, the energy 
absorption by the internal degrees of freedom is proportional to /^(T^). Another 
example is of quantum dots where parameters are varied adiabatically, like in 
pumping experiments, but with dots that are closed, so that their spectrum is 
discrete U. 

In Section g a specific RMT model is defined. For this model the ensemble 
average of the integral of the correlation function ( jl^ ) is calculated analytically. 
It is demonstrated that most of the contribution for long times originates from 
the nearest neighbor levels. In Section ^ the integral of the correlation function 
([12I), predicted by the nearest neighbor level spacing distribution, is calculated 
for various distributions, including some that are not related to RMT models. 
In Section ^ the results of this work are analyzed and discussed. 

2 Random matrix models 

The main purpose of this paper is to study ( [ill ) ^^"^ i^^ finite time integral (O) 
in the framework of RMT. The reason for this is that random matrices describe 
many characteristic properties of realistic quantum- chaotic systems [17, 25-28]. 



For simplicity our random matrices will depend on one external parameter X, 
so that we shall study the one dimensional version of dO), where the vector of 
parameters (R) is replaced by the scalar X. For each random matrix we shall 
be able to calculate both C{t) and I{t). 

2.1 A simple RMT model 

We wish to construct a random matrix model for some of the levels of a system 
whose quantum Hamiltonian depends on some parameter. The N levels we 
wish to simulate by the random matrix lie within an energy strip of width 
6E{N), that depends on N. Later on we shall be interested in studying the 
semiclassical limit. The meaning of taking this limit in the present context 
is to increase the density of levels in the SE-stiip: in the classical limit the 
spectrum becomes continuous. We shall work with the well studied Gaussian 
Ensembles pi| , |22|] . These are defined through four parameters: f3 which defines 
the symmetry of the random matrices, their dimension N, the mean value of 
their elements and their variance (given through the parameter /.P). All of these 
need to be chosen carefully in terms of parameters of the physical system, being 
simulated by the random matrix. The symmetry of the ensemble should be 
chosen to correspond to the real system. If the latter exhibits time reversal 
symmetry then the ensemble is the orthogonal one {(3=1). If the system does 
not exhibit this symmetry then the ensemble is unitary (/? = 2). The mean 
value of the matrix elements can be chosen to be zero, which corresponds to 
setting the ensemble average of the reference level (En) = 0. The mean level 
density satisfies the semi-circle law [|9[ ^ ^ : 

^-^^^1 otherwise ^ ' 



for large N, where we have used the definition: x = E/ ^APii^N. In subsection 
^.3| the relation between the parameters of the RMT model and the ones of the 
physical system will be discussed. 

We shall use the Hamiltonian introduced by Austin and Wilkinson |g] and 
model a parameter dependent system by the N x N random matrix: 

H{X) = Hi cos X + H2 sin X, (14) 

where Hi,2 are N x N random matrices from the same GOE or GUE ensemble. 
There are three advantages to working with H{X): (a) it belongs to the same 
ensemble that Hi 2 belong to, (b) the derivatives of its matrix elements belong 
to the same ensemble too because: 

dH{X)/dX = -HismX + H2COSX, (15) 

(c) the matrices H{X) and dH{X)/dX are statistically independent. If we insert 



H{X) and dH{X)/dX into Eq. 11, and then perform the ensemble average, we 
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Figure 1: The integral of the correlation function for GOE. Numerical results 
for TV = 3 (o), iV = 13 (v), iV = 53 (□) and N = 103 (a) are shown. Also 
shown is the large N approximation (Eq. |26| ) (line). The inset shows the long 
time behavior on a log-log scale. The number of ensemble members used is 10^ 
and the errors are of the order of A/i « 0.01. 




Figure 2: The integral of the correlation function for GUE. Numerical results for 
Af = 3 (o), iV = 13 (v), iV = 53 (V) and N = 103 (a) are shown. Also shown is 
the large N approximation (Eq. |29| ) (line). The inset emphasizes the long time 
behavior. The number of ensemble members used is 10^ and the errors are of 
the order of A/2 w 0.01. 



obtain: 



Cp{t) 



J2 \{dH{X)/dX\^ 



— {En — E„ 



(16) 



denotes RMT ensemble 



where {dH{X)/dX)^„^ = (n dH(X)/dX mj and , 

averaging. The subscript /3 will denote the symmetry: /3 = 1 for GOE and 
/3 = 2 for GUE. The correlation function Cp and its finite time integral Ip 
are ensemble averaged. The (. . .) will be dropped from these quantities for 
notational simplicity. 

The statistical independence of dH{X)/dX and H{X) implies 



Cp{t)= J2 (\{dHiX)/dXX 



— {En — En 



(17) 



while the fact that dH{X)/dX belongs to the same ensemble as H{X) implies 



(di/(X)/dX)„,„ 



(^(^))„,. 



/3m' 



(18) 



for m ^ n, leading to: 






— {En — Em) 



(19) 



We would like to make the connection between Ci3{t)/ (ifj? and the form 
factor: 



m-j 



1 



r{E) 



p{E + e/2p)p{E - e/2p) 



1 



e'^'^'^^de, 



(20) 



where p{E) — J^i ^{^i ~ E) is the density of states and p{E) is the smoothed 
density of states. The variable e is the energy measured in units of the mean 
level spacing l/'p{E) and r = t/Tu is time in units of the Heisenberg time, 
Th = h-p{E). 

Eq. hS can be written in the following form: 



Cp{T)/Pp^ = 



\^Ap{E + e/2-p)p{E - e/2-p)) - 5{e) 



-p\E) 



'de, (21) 



where the 5{e) results from the omission of the term to = n in the sum (|19|). 
Comparing the last equation with ( pO|) one can see that: 

(22) 



Cp{t)/Pp'=K{t)+5(t)~1. 
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In this work we are mainly interested in the time integral of the correlation 
function (fl): 



Ip{T)/pt?TH 



dT'Cf,{T')/(i^J.^ 



1/2- 



r dT'(^i-K{T') 



(23) 



In the limit r — > oo the term in the square brackets is just R^if- — 0), the two 
point spectral correlation function at zero energy separation. It vanishes as a 
result of level repulsion. 

In order to perform actual calculations we make use of the well known form 
factor for GOE and GUE [0 . It is standard to define 



6(t) := 1-K{t). 

For GOE it is given for example in Mehta's book (see |2^ p. 137): 

' 1 - 2t + r In [1 + 2r] t < 1 
-1 + Tln[i?±i] r>l ' 



6(r) 



from which one obtains: 
h{r) i 1- 



f^^TH 









-i)ln[l + 2T]] 



2r + l 
2r-l 



T< 1 
T> 1 



For T 



oo /i (r) falls off asymptotically as 

hir) 1 

fi^Tn 12t' 



For GUE (see [|| p. 95): 

6(r) = 

from which one obtains: 

hjr) j 



1-T 





T < 1 
T> 1 



r < 1 
r > 1 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



In order to compare the analytical results that hold in the infinite N limit 
with results for finite N, ensembles of the N x N matrices Hi and H2 of (|lj), 
belonging to GOE or GUE were generated numerically. The integral of the 
correlation function Ip was then calculated numerically, by ensemble averaging. 
The results are presented in Figs. |l] & || for GOE and GUE respectively and 
compared with ( |2^ ) and (p9|). Units where j3^^Th = 1 were used (see subsection 
.3|). The numerical errors in the figures were calculated according to 



A/^ 



iV,,, 



Ub 
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0.04 



Figure 3: Testing the nearest neighbor approximation (^) for GOE. Full nu- 
merical results for N = 13, including all level spacings (v), using only nearest 
and second nearest neighbors (x) and using only nearest neighbors (+). Also 
shown is the long time approximation, where only nearest neighbor spacings are 



taken into account (Eq 

10^ and the errors are of the order of A/i 



39) (line). The number of ensemble members used was 
0.005. 



where Nens is the number of matrices used in the ensemble average. 

For completeness the results for GSE are obtained with the help of (see [p2| 
p. 166): 



b{r) 



1 - ir + ir In I 1 - r I t < 2 



T> 2 ■ 

Following the calculation performed for the other ensembles one finds: 

h{r) r l^-[±{Ur-5r^) + l{r^-l)ln\l-T\] r<2 



4/z2rH 







r > 2 



(30) 



(31) 



2.2 Nearest neighbor spacing dominance and the long time 
limit 

The model (Q) is very specific in its dependence on the parameter X. An im- 
portant property of this model is the statistical independence between H{X) 
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0.08 - r 



0- ^f^ y »» ¥¥»ft^f»^ y y¥¥fft*»»^¥W¥¥» $$$ ^ $ »tH||$^ 



Figure 4: Same as Fig. pi but for GUE, compared with the long time approx- 
imation (|40| ) (hne). The number of ensemble members used was 10^ and the 
errors are of the order of A/2 ~ 0.003. 
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and dH(X)/dX . Such independence holds to a good approximation for disor- 
dered systems BG]. It is reasonable to make this approximation also for RMT 
models of chaotic systems. The reason is that most eigenstates look random, 
are statistically independent of the eigenvalues and therefore for many types 
of parametric dependencies the matrix elements of dH{X)/dX will look ran- 
dom and independent of the spectrum. Although this argument is reasonable 
for many types of parametric dependencies it is clearly not general. For the 
asymptotic behavior much less is required, since the long time asymptotics is 
dominated by the nearest neighboring levels. The reason for this dominance 
is that if T ^ 1 the terms in the sum (16) oscillate wildly as a function of 
m, so that the important net contribution is from the terms nearest to being 
stationary. These are obviously m = n± I. The approximation is therefore: 



C^i) 



{dH{X)/dX\ 



cos 



S(^" 



En-l) 



ioTt:^TH. (32) 



Only ?7i = n — 1 is required since we know that the matrix elements and eigen- 
value distributions are symmetric with respect to reflection around the middle 

2 

{dH{X)/dX)^_^^_^_^ by its mean value 
and ignore the contribution from its fluctuations resulting in: 



eigenvalue n. We further approximate 



Cp(t) « 2/3^^' 



cos 



■T {En — En-l) 

n 



(33) 



Hence we ignored 



AC^(t) 



{dH{X)/dX)_ 



f3fi' 



t 



{En ~ En-l) 



(34) 



in the t ^ Th limit. In the framework of RMT ( |33|) takes the form: 

/•OO 

Ci3{t)I(3h^ ^2 I ds Pp{s) cos (27rTs), (35) 



2 / ds Pp{s) cos {2'KTs), 
Jo 



where s is the nearest neighbor spacing in units of the mean level spacing A_E = 
l/p(0), the time r is measured in units of the Hcisenbcrg time Th and Pf}{s) is 
the probability distribution of s. The integral of the correlation function is: 



-^4w- ~ 2 r dr' Hds Pp{s) cos(27rr's) == - 
Pt^^Tn Jo Jo T^ 



ds -^^sm{2TrTs). (36) 



For the nearest neighbor level spacing distribution we use the Wigner surmise 
[ p!9| , that is exact for 2 x 2 matrices, and takes the form (see Eq. 202 in pi|| ): 



Pi{s) = 2*^^P 



(37) 
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for GOE, and 



32 

P2{s) ~ — 2S^exp 



-— s 

TT 



(38) 



for GUE. The integral ( p6| ) can be calculated for these distributions. For GOE 
one finds: 

ii|^ = 1 exp [~4ttt^] erf [z 2^t\ /z, (39) 

while for GUE one finds: 






2r exp 



(40) 



A crucial approximation in the long time regime is (|32D, where only the 
contribution of the nearest neighbors is taken into account. This approximation 
is tested in Figs. || & § for the model (|^. In the numerical test the nearest 
and next-nearest neighbor spacings are taken from the center of the matrix. We 
see from these figures that the approximation ( P4 ) is quite reasonable and it 
improves as time increases. The reason is that as time grows the oscillations of 
the various terms in (ITq) with energy become stronger, enhancing the dominance 
of the nearest neighbor contributions. The approximation is better for GOE 
than for GUE because the small s weight in the integral (^) is larger. 

The asymptotic behavior of erf \iy\li for large real y is exp[2/^]/(Y^j/) | |31[ . 
Therefore, we find for GOE: 



h(r) 



1 

47rT 



for T > 1, 



(41) 



and there is no characteristic time for the crossover from classical to quantum 
behavior. This decay is extremely close to ( p7| ) hence for large r the contribution 
of the nearest neighbor spacings accounts for nearly all of I\ (r) . 

For GUE the falloff is much faster since in the large t limit the Gaussian 
in (40) dominates, and there is a characteristic time 2/7r'^/^ (in units of the 
Heisenberg time), for the crossover to quantum behavior. For long time I^ij) 
is extremely small, while it is exactly zero according to (E9[) . Therefore also for 
GUE (p3) is an excellent approximation for the long time limit. 

For completeness let us present the results for the Gaussian symplectic en- 
semble (GSE) (for which /3 = 4), analogous to (|39| ) and (pl)|). For this we need 
to substitute the appropriate P/iis) into (pq). According to Eq. 202 in |2l|| : 



P,{s) 



)18 



367r3 



s exp 



64 



leading to: 



h{r) 



2 1 



3^ 
32 



r exp 



"64" 



(42) 



(43) 
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One can immediately see that there is a characteristic time scale for the quantum 
to classical crossover, as there was in the GUE case. As for the GUE case liir) 
is extremely small for long times, while (pw gives zero. 

In Figs. H & ^ it is demonstrated that the long time results of the simulation 
tend to the analytical formulas as more and more members are included in the 
ensemble. 

2.3 The short time Hmit and determination of the param- 
eters of the RMT model 

In order to make contact with a physical system one has to relate the RMT 
parameters /i and N with h and the parameters of the physical system. This 
is done in the short time limit. This limit is not universal and the dynamics 
of the chaotic system is not described by RMT. It is used only to determine 
the relation between the parameters. It will be assumed for concreteness that 
the system we wish to model by a random matrix is a two dimensional chaotic 
billiard (a free particle of mass m in a two dimensional box) of area A. The 
results of the paper do not depend on this assumption. First we establish a 
relation between the mean density of states of this model to the one of RMT. In 
the framework of RMT the semicircle law ( |l^ ) can be used for the mean density 
of states. If the density in the center of the strip of energies, modeled by the 
random matrix, coincides with that of the two dimensional billiard, 

TAo) = — = Vm^^7^2d{E)^Vm^^^- (44) 

The existence of the semiclassical limit for the correlation function (^Tj) [Q 
leads to another constraint. This constraint enables the expression fi in terms 
of h. It turns out that in the semiclassical limit the number of levels in a given 
interval grows with N. We shall give the explicit connection between classical 
parameters of the system, h, /i and N in what follows. For the model (Uj) and 
other models where statistical independence between dH{X)/dX and H{X) 
holds, (O) can be used. In this framework it is easy to take the semiclassical 
limit. In this limit the spectrum can be considered continuous for fixed time, so 
that the sum can be replaced by a suitably weighted integral: 

Cp{t) = Cp{t)~f3fi' (45) 

where 

C(3{t)^Pt? dxp^ix) cos -x/APfi^Nx ioTh^O, (46) 

Jo L" J 

and PfM^ is the contribution of the ra — n term in (|l9|). In the spirit of the 
RMT modeling we choose the eigenvalue En to be in the middle of the region 
described by RMT, therefore we set En — 0. The integral in (|4q) is known ||32|, 
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and we obtain: 

^ ^ for?i-^0. (47) 



C(3{t)^y/l3ij^-3i 



2^l3y?N- 



The correlation function has a characteristic time scale Tc — Ti / ^Jl3jJtFN . Since 
the integral (|l^) for //3(t) is convergent for all values of t ^ T^, it is well 
approximated by its value at < ^ Tc- Therefore, 7/3 (i — oo) = / is expected to 
take a classical value in the limit N ^ oo. One finds: 



/ = ny%iFN. (48) 

Now units where Am/2 = 2 and / = 1/2 are introduced. In such units fi and h 
are dimcnsionless and arc given in terms of N as: 

Pt^' = \n-^I\ (49) 

n = N-^l^. (50) 

The Heisenberg time in these units is: 

Th = /i -p^^{E = 0) - 4iVi/3, (51) 

and 

f^ii^Tu = 1, (52) 

while the characteristic time scale for the saturation of Ip to its classical value 
is Tc = 27V"2/3^ Finally, in these units: 

C/3(t) ~ ^^i^^^ ^ foriV->oo, (53) 

where r = t/Tn is the dimcnsionless time. The integral over the correlation 
function is. in these units: 

Ip{T)^ 1/2 -T, (54) 

and the approximation holds for r ^ 1. Now we can justify some of the assump- 
tions that we made. First of all, we see that the limit N ^ oo indeed corresponds 
to the limit h ^ 0. Secondly, the mean level spacing AE ~ 1/P2d(0) ^^"^ ^c 
decay to zero (as N^^^^) in the limit N -^ oo, as expected. 

3 The long time behavior predicted from the 
nearest neighbor level spacing distribution 

In the previous section the integral of the correlation function Ip (r) was studied 
for RMT models. One conclusion was that it is dominated by the nearest 
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neighbor level spacings. It was found also that there is a big difference between 
GOE on the one hand and GUE and GSE on the other. In this section we shall 
study the decay of the correlation function if the nearest neighbor level spacing 
distribution is given and will not rely on the assumption of an invariant RMT 

2 



[dH{X)/dX)^^^_ 



are 



ensemble. We will also assume that the fluctuations of 

not important and this quantity can be replaced by the constant /3/i^. It will 
be assumed that the distribution of the nearest neighbor level spacings is of the 
form: 

P/3(s) =cs^exp[-as2] (55) 

where a and c are constants. The RMT distributions (p7|), ( {38|) an d (42| ) for 
(3 — 1, 2, 4 are of this form. The integral for the correlation function (^) takes 
the form: 



Ip{t)/P^?Th ^- ds s^^^exp [-as'] sin(2^rs) = — --J^(y) 



where 



/•OO 

Xfj{y) = I ds s'^~^ exp [— s^] sinsy, 
Jo 



(56) 



(57) 



with y = 2-KT j \fa. The decay of this function will be explored in what fol- 
lows, for arbitrary /3. One can verify that this function satisfies the ordinary 
differential equation: 



dy- 



,Ip{y)+y^i-Ip{y) + ^My)^0 



1^0) = 0; Ax^(0) = ir 



13+1 
2 



Making the substitution: 



I/3(y) = //3(2//\/2)cxp[-yV8] 



(58) 



(59) 



and changing to the variable x — j//v2, one arrives at a new differential equa- 
tion: 

r /"(x)-[xV4+l/2-/3]/(x)=0 



p+i 

2 



(60) 



\ /MO) = 0; /^o) = ir 

Eq. pOl is a well known equation, and its solutions are Parabolic Cylinder Func- 
tions^: U [1/2 - /3, x] , V [1/2 - fi.x]. For arbitrary /3, the solution of (|o|) is 
a linear combination of these two functions: 



fp{x)^Ap U 



'1 

2 


-/3,x 


+ BpY 


'I 
2 


-/3,x 



(61) 



in which Ap^Bp arc constants to be determined from the initial conditions. In 
particular, if (3 is an odd natural number it turns out that Ap — 0, if it is an 
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even natural number B/j — while for non-integer /3 both Ajs and Bp are non- 
vanishing. This fact is very important for the asymptotic behavior oiX/jly), 
which is determined by the large x behavior of U [^ ~ P,x\ ,V [^ — (3, x\ pl[: 



2-/3,0. 



r/3-l 



(3,x 



exp [-a;^/4] 



(62) 
(63) 



as a; ^ cxD. One can immediately deduce that the first solution is subdominant 
for all /3 except for the special case when it is an even natural number, for which 
Bf} = in (61). Since 1j3{y) is proportional to Ii3{t) and y is proportional to r, 
for any fixed (3 ^ 2n (n = 1, 2, 3, . . .): 



as r — > oo, while for P — 2n (n — 1, 2, 3, 

7/3 (r) ~exp 



-l3 



./5-1 



(64) 



(65) 



as T — > oo. This is precisely the type of behavior found for the random matrix 
ensembles treated explicitly (41, ^ & p3| ). 

Finally, one wonders what would the analogous results be in the case of the 



Poisson distribution. Returning to (32), one obtains for the correlation function: 

POO 

CpiT)Ki2ap / (is Pp(s) cos (27rTs), 
Jo 



(66) 



where ap is the variance of the off-diagonal matrix elements between nearest 
neighbor levels, which we shall leave unspecified, as we are only interested in 
the behavior as a function of r. An approximation similar to the one leading to 
( p3| ) was made. For the Poisson case the nearest neighbor spacing distribution 
is: 

Pp{s) = exp[-s], 

leading to: 

2a% 



Cp{t) 



1 + {2t:tY' 



(67) 
(68) 



The integral over this expression is: 



Ip{r) 



(jpTn 



arctan (27rT), 



the asymptotic behavior of which is given by: 



Ipir) 



1 

2ttt 



ctITh/2, 



(69) 



(70) 
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as r -^ oo. In the absence of level repulsion one indeed finds that the integral of 
the correlation function does not vanish. This does not result in any contradic- 
tion with the classical limit where for an integrable system I(t = oo) = 0. The 
reason is that for integrable systems the eigenfunctions of neighboring energy 
levels typically have an exponentially small overlap (in l/Ti). This small overlap 
is the physical reason for the Poisson distribution. Therefore in the classical 
limit CTp — > exponentially fast in l/h. 

Many systems that are neither integrable nor chaotic were found to have a 
semi-Poisson distribution |Q. For this case, that shows linear level repulsion, 
the nearest neighbor spacing distribution is: 

Psp(s) = 4sexp[-2s], (71) 

The integral of the correlation function was calculated along the lines of the 
calculation for the Poisson distribution. The result is: 

/5P(r)«2a^pT«^^^. (72) 

It decays like l/r in the long time limit. For the semi-Poisson distribution there 
is level repulsion and indeed the integral of the correlation function decays with 
time. 



4 Summary and Discussion 

The correlation function of the force applied by a fast quantum system on a 
slow classical one is calculated within the leading order correction to the adi- 
abatic approximation following Berry-Robbins and Jarzynski. Its finite time 
integral I{t) of (|l^ ) is proportional to the dissipation rate on the time scale t. 
In the present work /(t) was calculated under various statistical assumptions. 
In Section g it was studied in the framework of RMT. For the specific depen- 
dence of the model ( [l4|) on the parameter X, the Hamiltonian H and dH/dX 
are statistically independent [gj. For this case it was shown that up to a pro- 
portionality constant, the integral of the correlation function is simply related 
to the integral of the form factor (by Eq. ^3|) . Since the form factor is known 
in RMT, the integral of the correlation function was calculated and found to 
vanish for times beyond the Heisenberg time for GUE and to fall off as a power 
law for GOE. This is a remarkable and surprising difference. The result for GSE 
is similar to the one found for GUE, except that the integral of the correlation 
function vanishes after twice the Heisenberg time. The properties of the model 
( |l4|) are satisfied approximately by many systems [pO[ , therefore it is expected 
that the results of this work are relevant for a wide range of problems. For the 
model (O) we have shown that for long times the results are dominated by the 
nearest neighbor spacings. If only these spacings are taken into account one 
finds that for long times the integral falls off as a power law for GOE and as a 
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Figure 5: Testing the dependence on the size of the ensemble for GOE with 
A^ = 13 where aU level spacings are taken into account: 10'' members (v), 10^ 
members (*) and 10® members (+). The errors are of order A/i « 0.01, A/i « 
0.005 and A/i « 0.001 respectively. Also shown is the long t approximation 
(Eq. M) (hne) 
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Figure 6: Testing the dependence on the size of the ensemble for GUE with 
A^ = 13 where aU level spacings are taken into account: 10"* members (v), 10^ 
members (*) and 10^ members (+). The errors are of order A/2 w 10^^, A/2 ~ 
lO"** and A/2 ~ 10~^ respectively. Also shown is the long t approximation 
(Eq. ^) (solid line). The errors are marked by dotted lines, the top line is for 
10"^, the middle for 10^ and the bottom for 10^. 
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Gaussian, where the characteristic time is proportional to the Heisenberg time 
for GUE and GSE. These results do not require all the properties of the model 
(Oh. They require only that the contribution of the fluctuations in the absolute 
value of the matrix elements between nearest neighboring states (neighboring 
in energy) is negligible, namely that the contribution of AC^(t) of ( p4[ ) can be 
ignored compared to the one of C^(i). This is clearly a weaker assumption than 
complete statistical independence between H and dH/dX. Therefore we expect 
the difference found between GOE and GUE (as well as GSE) to be generic for 
RMT models with various dependencies on the parameter X. The long time 
behavior of the RMT models is expected to provide a faithful representation of 
the behavior of chaotic systems, since it is dominated by the small level spac- 
ings. For short times, on the other hand, the behavior depends on the specific 
properties of each system. The short time behavior of the RMT model was pre- 
sented here only to set the relation between the constants of the RMT model 
and the ones of the chaotic system. 

The assumption of the dominance of the contribution of nearest neighbor 
level spacings, together with the assumption that ACff{t) of Q) is negligible 
enables the calculation of the integral I{t) of (|l2|) for various distributions of 
nearest neighbor level spacings even if these do not necessarily originate from 
RMT models. For the distribution (pa), that is a generalization of the distribu- 
tions found for GOE, GUE and GSE, one can calculate //^(i) for various values 
of p. In Section S it is found to decay as t~^ for all values of /3, except when /3 is 
a positive even integer, for which it decays like a Gaussian with a characteristic 
time that is proportional to the Heisenberg time. What is special when /3 is a 
positive even integer? For these values the integral ( |56| ) can be extended to the 
range [— oo, cx)]. The integrand is an entire function, the contour of integration 
can be deformed in the complex plane and the integral is dominated by a saddle 
point. For other values of f3 an extension of the integral to negative t, so that 
the integrand is analytic, is impossible. The point s = is an end point and 
for large t the integral is dominated by it, leading to power law decay. For 
non- integral f3, the point s = is also a singular point. It would be nice to find 
a more physical explanation for this difference between the various ensembles. 
For completeness the integral of the correlation function was calculated for the 
Poisson and the semi-Poisson distributions. 

The various RMT formulas (Eqs. |2g & ^9|) are developed for the limit of 
infinite matrices. This limit is approached extremely fast, as can be seen in 
Figs. & H The convergence to the average, as a function of N^ns, the number 
of members of the ensemble, is slow (see Figs. 1^ & |6|). 

The crucial approximation that was made generalizing the results beyond the 
model ( [l4|) was neglecting ACfj{t) of (|34|). Although reasonable, its validity for 
chaotic systems should be checked. The results may hold also for mixed systems 
if sticking to integrable regions does not take place on time scales relevant for 
the calculation. For chaotic systems corrections of order higher than the leading 
one, in the adiabatic approximation, may lead to different behavior after some 
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time {T2 or Tlz)- Dephasing, as a result of the coupling to the environment 
will destroy the quantum correlations on a time scale T^. In experiments of 
the type mentioned in the end of the Introduction, dephasing is always present, 
and if T^ ^ T2,Tlz, the system will dissipate energy at a rate proportional to 
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